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1 Introduction 



Most of the concepts used to understand the viscoelastic properties of chemical 
and physical gels of flexible polymers require the persistence length ip to be 
significantly smaller than other characteristic scales such as the filament length L, 
the distance between crosslinks or the width of reptation tubes. This condition no 
longer holds for networks of semiflexible polymers. One prominent family of such 
polymers are cytoskeletal biopolymers like F-actin, intermediate filaments and mi- 
crotubules. An impression of the typical conformations of the filaments and the 
relative magnitude of the characteristic length scales can be gained upon inspec- 
tion of the following electron micrograph of a semidilute actin solution. The most 
striking features of these networks are the enormous length and relatively elongated 
structures of the constituent biopolymers. Actin filaments have a diameter of 7nm 
and can reach lengths up to 30-100 fim in vitro 0, and several microns in vivo. 
The persistence length is approximately 17yum [H, H, P], quite large compared to 
typical distances between neighboring filaments which is in the range of a few tenth 
of a micron. This combination of length scales allows biopolymers to form networks 
at very low volume fraction, so that solutions of less than 0.1% volume fraction of 
polymer are still strongly entangled. Thus only a small amount of material needs to 
be produced by the cell in order to generate a sufficiently strong network. This fact 
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Figure 1: Left: Electron micrograph of a 0.4mg/ml actin solution. The bar indicates 1 
/im. Right: Sketch of the wormlike chain as a space curve r(s). 



is not only of considerable biological relevance but also facilitates interpretation of 
dynamic light scattering experiments and their relation to theory |0, H. 

There are several motivations f^, [1^, |Tl], |I2| for studying the viscoelasticity of 



cytoskeletal networks H, |T|, [T5|, |16|, |T7|, |T|, l^, |2T|, |2|: 

(I) From a polymer physics perspective they are interesting because their be- 
havior is expected to be determined by principles and mechanisms different from 
those established for flexible networks. In fact much of the physics behind the 
viscoelasticity of semiflexible polymer networks is only being explored recently 

nn, m, ii,i3,|8|,|9ii3. 

(II) From an experimental point of view, semiflexible polymer networks are inter- 
esting because several semiflexible polymers have persistence lengths on the order 
of several /xm or even mm. Thus techniques such as optical microscopy of single 
fluorescence labeled filaments or attached colloidal probes can be used to study the 
behavior of the network at the single polymer level [|1^ 
polymers this has been possible only in simulations (see, e.g.. 



19, 22, 311. For flexible 



(III) Finally, the viscoelastic properties and regulation of semiflexible polymer 
networks both inside cells and in the extracellular matrix are of signiflcant im- 
portance for the mechanical stability and properties of biological tissue, for cell 
locomotion, adhesion and force generation |1^, pJ]] . 



The Cytoskeleton: Structure and Biological Role 

Living cells need both the ability to maintain their shape when exposed to shear 
stresses exerted by their active contractile machinery or by fluid flow in blood ves- 
sels and the ability to reorganize their shape and internal architecture as is the case 
in cell migration and mitosis. The structure responsible for the mechanical and 
dynamic properties of the cell is the cytoskeleton, a rigid yet flexible and dynamic 
network of proteins of varying length and stiffness. Most eukaryotic cells contain 
three types of protein fllaments comprised of actin, tubulin and intermediate filament 
proteins such as vimentin. These, as well as the plasma-membrane associated flla- 
ments make up the cytoskeleton |^ . There is also a range of accessory proteins for 
each of the cytoskeletal fllaments which allow for control of nearly all mechanically 
relevant properties of the protein fllaments [^, |3^. Let us just mention one exam- 
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pie, gelsolin, which is frequently used in rheological experiments. It caps the end of a 
growing actin filament and can thus be used to regulate the average filament length 
in actin solutions. There are many other proteins with different tasks ranging from 
initiating and terminating polymerization over introducing crosslinks and forming 
lateral arrays of filaments to even changing the stiffness of the filaments. 

An example of the biological role of the cytoskeleton is the migration of an 
amoeba. Its motion is initiated by adhesion driven spreading of the cell membrane 
on the substrate followed by gelation of actin in the advancing lobe (pseudopodium). 
The cycle is completed by retraction of the rear end and a gel-sol transition or fiber 



formation in the advancing front [Tl|. This process is of course quite complex and 
there is a subtle interplay between regulatory mechanisms and the material proper- 
ties of the cytoskeleton. But, understanding the basic physical principles determin- 
ing the viscoelasticity of actin networks is certainly a prerequisite in understanding 
such a biological process. In the following we will address the following questions: 
(i) Can we even understand the physics of a one-component system such as a pu- 
rified F-actin solution (depicted in Fig. P? (ii) Can we identify the basic physical 
principles underlying the observed viscoelastic behavior? (iii) How is it different 
from the physics of long flexible coils? 



2 Single Chain Properties 

The model usually adopted for a theoretical description of semiflexible chains is 
the wormlike chain model [|l], Here one describes the filament as a smooth 

inextensible line r(s) of length L parameterized in terms of the arc length s. The 
statistical properties are determined by an effective free energy (the "Hamiltonian" ) 



mr{s)}) = ^ f^ds 



which measures the total elastic energy of a particular conformation by the integral 
over the square of the local curvature weighted by the bending modulus k. The 
inextensibilty of the chain is expressed by the local constraint, |t(s)| = 1, on the 
tangent vector t(s) = dr/ds. We will see that this constraint is essential for a correct 
description of the static as well as the dynamic properties of semiflexible polymers. 
Due to the mathematical complications resulting from the inextensibilty only few 
of the statistical properties of the wormlike chain can be extracted analytically, the 
best known being the exponential decay of the tangent-tangent correlation function 
(t(s)t(s')) = exp (— |s — s'\/ip) with the persistence length £p = n/kBT, the mean- 
square end-to-end distance [|| TZ^ := {[r{L) — r(0)]^) = LF' f o{L / C-p) , and the radius 
of gyration ^3 7^2 = £2 {f^{L/£p) - 1 + L/3£p), where /d(x) := 2(6"^ - 1 + x)/x 



2 



is the Debye function. 



2.1 Force- Extension Relation 



One of the most obvious differences between flexible and semiflexible polymers is 
their response to external forces (see Fig. ^ (left)). In the flexible case the response 



Physics of Solutions and Networks of Semiflexible Macroniolecules 



4 



is isotropic and proportional to l/ksT, i.e., the Hookian force coefficient is pro- 
portional to the temperature. When the persistence length is of the same order of 
magnitude as the contour length, the response becomes increasingly anisotropic. 



X d = 3, 40 segments, L/t^ = 10 0.5 



//m 

entropy 
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Figure 2: Left: A flexible chain's linear response is purely entropic and isotropic. The 
elastic response of a stiff rod is extremely anisotropic. Transverse forces (bending) lead 
to a purely mechanical response, whereas the longitudinal response (buckling) is charac- 
terized by a spring constant inversely proportional to the temperature. Right: Numerical 
results for the end-to-end distribution function of a (discretized) wormlike chain in d = 3 
dimensional space (numerical data from Ref. [^8|). With increasing stiffness there is a 
pronounced crossover from a Gaussian shape to a form with the weight shifting towards 



full stretching. The dashed line indicates the Daniels approximation [39|. 



Then the linear response of the chain depends on the orientation of the force 
with respect to the tangent vector at the clamped end. Transverse forces give rise to 
ordinary mechanical bending of the filaments and the transverse spring coefficient is 
proportional to k. The linear response for longitudinal forces is due to the presence of 
thermal undulations, which tilt parts of the polymer contour with respect to the force 
direction. The effective longitudinal spring coefficient turns out to be proportional 
to k'^/T indicating the breakdown of linear response at low temperatures (T 0) or 
very stiff ffiaments {ip — > oo). This is a consequence of the Euler buckling instability. 
Note also that for the special boundary conditions of a grafted chain (as depicted 
in Fig. (left)) the linear response of the chain can even be worked out exactly 



for arbitrary stiffness p^; these calculations use the fact that the conformational 



statistics of the wormlike chain is equivalent to the diffusion on the unit sphere p 



2.2 Radial Distribution Function 

An important quantity describing the statistical properties of the chain is the 
probability distribution of the end-to-end vector G{r;L) = (5(r — R)). For a 
freely jointed phantom chain this function is known exactly As for any 

model with short-ranged interactions it converges quickly to a Gaussian distribution 
GqIt; L) ~ exp (— 3r^/4£pL) for an increasing number of segments. For chains that 
are at least some 10 £p long the Gaussian can serve as an excellent approximation to 
G(r; L) for many purposes. For the freely jointed chain the persistence length is 
independent of temperature because its microscopic origin lies in steric constraints 
rather than in the bending stiffness of the backbone. 
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For a Gaussian chain, the separation by a given distance r of any two seg- 
ments with preferred mean-square distance 2ipS is punished by the free energy 
cost F[r) = — ksT In Goir; s) = const. + SksTr'^ /AipS quadratic in the end-to- 
end distance. Due to the Euler instabihty this is very different for semiflexible 
chains. The characteristic feature of the physics of beam buckling is that the en- 
ergy Eel of a straight rod is an almost linear function of its end-to-end distance 
R, Ed ~ fc ■ {L — R). Here fc = ktt^/L^ is the critical force for the onset 
of the Euler instability. Neglecting fluctuations around the classical contour this 
would lead to an end-to-end distribution function with maximum weight at R = L, 
G(r; L) oc exp[— /c ■ (L — r)/kBT]. Note that with such an approach we completely 
ignore entropic effects which are the only contributions in case of the freely jointed 
chain, discussed above. In order to correct for this omission we have to multiply 
the above Boltzmann weight by the relative number of allowed conformations. This 
becomes most obvious for a completely stretched chain, where up to global rotations 
only one possible configuration exists and consequently the end-to-end distribution 
function has to vanish. These qualitative arguments lead to the shape of the distribu- 
tion function shown in Fig. ^ (right). The actual form of the end-to-end distribution 
function can be obtained within a quantitative analysis of the wormlike chain. 



2.3 Dynamic Light Scattering 

A useful experimental technique for investigating the short time dynamics of semi- 
flexible polymers is dynamic light scattering (DLS). In DLS experiments one directly 
observes the dynamic structure factor (/(q, t). We focus on the ideal case of a dilute 
or semidilute solution of semiflexible polymers, where the scattering wavelength is 
much smaller than the mesh size. We also assume a separation of length scales, 
a <^ X < ip, L, i.e., the scattering wavelength A is large compared to the monomer 
size a but small compared to the characteristic mesoscopic scale defined by L and ip. 
As a consequence the contributions to the time decay of (^(q, t) from center of mass 
and rotational degrees of freedom of the chain are strongly suppressed as compared 
to contributions from bending undulations. Moreover, for this case [^, |^, ^] 

the structure factor can be written as exp(— g^r^(t)/4) with the local mean square 
displacement r'j_{t) ~ t^^^: 

^(q,t)ocexp[-(7,t)3/4], (2) 
where 7^ ~ Q^^"^ / C±(p^'^ [0) 0- Such a stretched exponential behavior has been con- 



firmed experimentally with very high accuracy for F-actin solutions |^^. However, 
a more careful analysis reveals that it cannot hold for very short times. For times 
shorter than C±/K,q'^ the bending forces can be considered weak and the contour 
obeys the fast wiggling motion imposed by hydrodynamic fluctuations. As a conse- 
quence the initial decay is of the form g{q,,t) oc exp(— 7g°'*t) with [0] 

For polymers, which are not quite as stiff as actin, e.g. for so called intermediate 
filaments, this initial decay regime is readily observed in light scattering experiments. 
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Figure 3: Left: The g— dependence of the initial decay rate compared to DLS data on 
actin 1 4^. Right: Comparison of the classical result for a swollen Zimm chain ||4^ with 
Eq. (|) and quasi-elastic light scattering experiments with the semiflexible biopolymer 
fibrin [^]. The data were most kindly provided by G. Arcovito (see also this volume). 



A very convincing confirmation of these theoretical results has recently been found 
in fibrin systems (see Fig. Analyzing the data by Eq. ^ allows one to 
estimate the friction coefficient (± entering the Langevin equation or, equivalently, 
the thickness a of these filaments. 



3 Collective Properties 

In conventional polymer systems made up of long flexible chain molecules the vis- 
coelastic response is entropic in origin over a wide range of frequencies |^6|. For 



semiflexible polymers a complete understanding of the viscoelastic response is com- 
plicated by several factors. First of all, there are several ways by which forces can 
be transmitted in a network. This can either happen by steric (or solvent-mediated) 
interactions between the filaments or by viscous couplings between the filaments and 
the solvent undergoing shear flow. It is a priori not at all obvious which if any of 
these coupling will dominate. In the case of flexible polymers it is generally believed 
that macroscopic stresses are transmitted in such a way that these transformations 
stay affine locally, i.e. that the end-to-end distance of a single fllament follows the 
macroscopic shear deformation [Q. Second, single fllaments are anisotropic elas- 
tic elements showing quite different response for forces perpendicular or parallel to 
its mean contour. Therefore one has to ask what kind of deformation of the actin 
filament is the dominant one and whether due to the anisotropy of the building 
blocks of the network macroscopically affine deformations stay affine locally. In the 
following we will address some of the issues raised. 



3.1 Plateau Modulus for Entangled Solutions 

If solutions of semiflexible polymers are sufficiently dense and are probed on suffi- 
ciently short time scales (typically in the range of 10~^ Hz to 1 Hz) they will exhibit 
a "rubber plateau". Its existence is in general traced back to a time scale separation 
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between the internal dynamics and the center of mass motion of the polymers. An 
externally imposed shear stress will then be transmitted to the individual strands, 
whose response will determine the magnitude of the modulus. This many chain 
problem is usually reduced to a single chain model by making certain assumptions 
on the effect of the mutual steric constraints on the conformation of a single filament. 

In what might be called the affine model the "phantom model" is adopted to 
semiflexible polymer systems [^. It is assumed that upon deforming the network 
macroscopically the path of a semiflexible polymer between two entanglement points 
is straightened out or shortened in an affine way with the sample. The macroscopic 
modulus is then calculated from the free energy cost associated with the resulting 
change in the end-to-end distance. Since in a solution forces between neighboring 
polymers can only be transmitted transverse to the polymer axis and there is no 
restoring force for sliding of one filament past another, it is however hard to imagine 
that entanglements are able to support longitudinal stresses in filaments. The mod- 
ulus predicted in the affine model should scale as oc c^^/^ and leads to absolute 
values of the order of 10 Pa; such high values are at odds with the low values ob- 



served in recent experiments on F-actin solutions ||4^. It was therefore argued [|T^] 
that such models are more appropriate for crosslinked networks, where they would 
predict a plateau value G° ~ kBTip/C,^. But, even in such chemical networks with 
crosslinks present it is a priori not obvious that local deformations on the scale of a 
single filament are actually affine and that longitudinal stresses in the filaments are 
the dominant contribution to the plateau modulus (see also section |3.3|) . 

Recent theoretical and experimental studies [^, ^ based on Refs. |]SD|, ^] 
suggest a different view. Here one considers the free energy cost of suppressed 
transverse fluctuations of the polymers that comes about by an affine deformation 
of the tube diameter. According to Odijk |^ the mean distance between collisions of 
a tagged polymer with its surrounding tube with diameter d is given by the deflection 
length Le Since each of these collisions reduces the conformation space 

there is a free energy of the order of kBT. The total free energy of u = c/L 
polymers per unit volume becomes F ~ v ksT L/ L^.. To be able to compare these 
results to experiments one needs to know how the tube diameter d depends on the 
concentration of the solution or equivalently on the mesh size := ^J'i|vL. In 
other words we have to determine the average thickness (i of a bend cylindrical tube 
in a random array of polymers as depicted in Fig. ^ (left). 

The contour and thickness of the tube will be determined by a competition 
between bending energy favoring a thin straight tube and entropy favoring a curved 
thick tube. This competing effects define a characteristic length scale which turns 
out to be Lg. For length scales below Lg the tube will be almost straight and we can 
estimate its thickness as follows. Upon restricting the orientations of the polymers 
to being parallel to the coordinate axes the density of intersection points (black 
dots in Fig. ^ (left) will be 1/^^. Hence for a tube of length Lg the line density of 
these intersection points projected to a line perpendicular to the tube increases as 
-^e/^m which implies that the tube diameter decreases with increasing tube length 
as (i ~ ^^/Lg. Hence one finds Le = (Cm^^^)^^^ leading to the following form of the 



Physics of Solutions and Networks of Semiflexible Macroniolecules 



8 




0.3 0.5 0.7 1 

c [mg/ml] 



Figure 4: Left: A semiflexible polymer can trade bending energy for a wider tube. The 
configuration of the constraining polymers (dots) is the same as in the upper figure. Right: 
Comparison of the predicted shear modulus ||5^ to experiment The total length L 
and persistence length £p were set to ip = 17 iJ,m and L = 16 //m, respectively pi]]. 



free energy and hence the plateau modulus 

G°~F~A;bT£;^/^c^/^ (4) 

The above scaling law is included as a limiting case in a more detailed analysis 
concerned with the calculation of the absolute value of the plateau modulus . 



The same scaling result has been obtained previously |^ using a different scaling 
argument. 

Recent experiments seem to favor the above tube picture, where the plateau 
modulus is thought to arise from free energy costs associated with deformed tubes 
due to macroscopic stresses. Fig. ^ (right) shows the results of a recent measurement 



of the concentration dependence of the plateau modulus in F-actin solutions 
which agrees well with the scaling prediction G° oc c^/^. 

3.2 Viscoelasticity and High Frequency Behavior 

At frequencies above the "rubber plateau" (i.e. above 1 Hz for a typical F-actin 
solution) a power-law increase of the storage and loss modulus with frequency, 
G'{uj) oc G"{uj) (X uj^l^ , has been observed |Q |55], |5^. It is tempting to speculate 



that it is somehow tightly connected with the anomalous subdiffusive behavior of the 
segment dynamics of a single filament. But in view of the actual micro- rheological 
experiments, where one observes the mean-square displacement of a bead of diame- 
ter larger than the mesh size and hence couples to a large number of filaments, it is 
not obvious how this comes about. A thorough understanding would need to explore 
the nature of the crossover from local dynamics dominated by filament undulations 
to the collective dynamics of the network and the solvent. 

At present there are two different theoretical approaches based on different as- 
sumptions on the nature of the dominant excitations of the individual filaments 
generated by the beads embedded in the network. In one class of theoretical models 
one takes over the above mentioned "phantom model" to the high frequency behav- 
ior |]57], p8|. It is assumed that under an applied shear deformation the filaments 
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undergo affine deformations on a length scale of order Lg implying longitudinal 



stresses on single filaments. In the high frequency regime this leads to [57, 5S 



G*{u) = ^u{kBTy/x/\zuja)-'/\ (5) 

independent of the entanglement length Lg. 

A complementary theoretical approach ||5^ starts from an effective medium de- 



scription for the polymer solution at large scales which crosses over to the single 
polymer picture at about the tube diameter. The low frequency response is due 
to peristaltic modes of the effective medium. At high frequencies, the penetration 
depth for these modes falls below the tube diameter and the excitations are bound 
to the polymer backbones. Assuming that the forces between polymers are trans- 
mitted by binary collisions, the transverse modes that make up the plateau modulus 
according to the tube model, are also responsible for the high frequency response. 
This again leads to an cu^/^— asymptotics of G*{uj) at high frequencies. However, the 
model describes the crossover to and the moduli within the plateau region and allows 
scaling predictions for the relationship between plateau modulus and entanglement 
frequency. The information contained in the viscoelastic moduli is conveniently 
expressed in terms of the density of relaxation modes. Preliminary investigations 
show that already the simplest scaling assumption for this density (which certainly 
greatly over-simplifies the complicated crossover from single polymer dynamics to 
the effective medium modes) leads to excellent agreement with experimental data. 
Using the fluctuation-dissipation theorem, the long time behavior of the dynamic 
structure factor in semidilute solutions can also be derived. 

Which one of these theoretical models is capturing the correct physics is not clear 
at present. It may well be that the actual physical mechanism is different from both. 
There is certainly a tremendous need for more detailed experimental studies which 
not only measure the power-law dependence of the modulus but also determine the 
concentration dependence of the entanglement frequency. 

3.3 Effect of crosslinking 

For a crosslinked network of semiflexible polymers bending and compressing forces 
can be transmitted to the filaments. Both for networks where the mesh size is 
very small compared to the persistence length so that the longitudinal elastic re- 
sponse of the polymers is dominated by their Young's modulus and for networks 
with larger mesh size where thermal undulations are crucial in understanding the 
elastic response of single filaments compression is a much stiffer mode of 

deformation than bending. Unless highly ordered network geometries are assumed, 
it is not clear which of the two modes will dominate the elastic response. Different 
assumptions on the real or effective network geometry can either favor the bend- 
ing modes as in or the compressional modes as in leading to substantially 
different predictions for the modulus. 

We used a two-dimensional toy model to investigate which type of deformation 
mode is dominant in a disordered crosslinked network. Sticks of length L were 
placed randomly on the plane and crosslinked at every intersection with another 
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stick. Crosslinks were inextensible. Sticks were assigned a Young's modulus E 
and a diameter r resulting in force constants fccomp = vrr^i? for compression and 
^bend = Sur'^E/SL'^ = 3k/L^ for bending the rod with one end clamped. Units were 
chosen such that L = 1 and n = 1. The model was subjected to periodic boundary 
conditions, strained and the linear elastic response calculated by the method of 
finite elements. While this is a purely mechanical model it captures the essential 
features of two very different force constants and disorder. Entropic contributions 
from fluctuations of the crosslink positions are not expected to be significant for 
dense networks. 

Which of the two modes dominates the elastic behavior was determined by keep- 
ing /cbend fixed and varying /ccomp- We observe that for slender rods or low densities a 





Figure 5: Left: Network of sticks for p = 50, L = 2 and a = 0.01. The color code indicates 
the load distribution with energy decreasing from red to blue. Right: Dependence of the 
shear modulus on the ratio A:comp/^T/^bcnd for networks with L = 15. 



certain point the modulus ceases to depend on fccomp, indicating that the elasticity is 
dominated by bending modes. While these two-dimensional results are certainly not 
straightforwardly applicable to three-dimensional networks we will nevertheless try 
to get a feeling for the scales involved. Network densities can be compared roughly 
by using the average distance Lc between intersections as a measure: A cytoskeletal 
network might have Lc ~ 0.1 fim with typical filament lengths of 2 fim correspond- 
ing to a two-dimensional density of p ~ 20 and an aspect ratio of a ~ 0.002 resp. 
^comp/^bend ~ 10^. Comparisou with Fig. | shows that this would just place the 
network in the bending dominated regime. This might, however, be different for 
different scales or if more order is present in the network than assumed here. For a 



more detailed analysis of the random stick model see Ref. 60 



4 Summary and Future Perspectives 

We have seen that the cytoskeleton is a composite biomaterial with a wide variety of 
interesting viscoelastic properties. In particular F-actin solutions and networks pro- 
vide a model system for a polymeric liquid composed of semiflexible polymers which 
is accessible to a complementary set of experimental techniques ranging from direct 
imaging techniques over dynamic light scattering to classical rheological methods. 



Physics of Solutions and Networks of Semiflexible Macroniolecules 



11 



From these studies it has become quite obvious that semiflexible polymer networks 
require new theoretical models different from conventional theories for rubber elas- 
ticity. The nature of the entanglement in solutions of filaments is very different 
from flexible coils. In a frequency window where an elastic plateau is observed a 
tube picture where the modulus results from the free energy costs associated with 
the tube deformations seems to be sufficient to explain the observed concentration 



dependence of the plateau modulus and even its absolute value ^3 



Outside the rubber plateau in the high-frequency as well as the low-frequency 
regime the situation is less clear. Micro-rheology and dynamic light scattering ex- 
periments allow us to access the short-time dynamics of the filaments within a net- 
work. Here a theoretical model which describes the combined dynamics of network 
and solvent in this regime is still lacking. At present there are two quite different 
approaches which either start from a continuum medium approximation or from 
a single- filament picture. Obviously both are just limiting cases and a molecular 
theory needs to explain how starting from the single-filament dynamics including 
interactions with the solvent and the neighboring filaments leads at some length 
and time scale to collective behavior, which might be described in terms of some 
continuum model. 

Another very important question is concerned with the effect of chemical 
crosslinks on the mechanical properties of semiflexible polymer networks. This is 
of prime interest for both cell biology and for polymer science. In cell biology one 
would like to know how the material properties (e.g. elastic modulus, time scales 
for structural rearrangement and stress propagation) change with the network ar- 
chitecture and the mechanical and dynamic properties of the crosslinks. From the 
perspective of polymer science it connects cytoskeletal elasticity with the very active 
fields of transport in random media and elastic percolation. In section we have 
presented a numerical study using a two-dimensional toy model. One can certainly 
not expect that such a simplified model leads to quantitative results, but we think 
that some of its main features carry over to the more complicated situation of a 
three-dimensional network. 
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